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O ' ABSTRACT 

co 

qq , Aims. Our aim is to study the production mechanism for very-high-energy (VHE; >100GeV) y-rays in distant active galactic nuclei (AGN) and 
»j . use the observed VHE spectrum to derive limits on the Extragalactic Background Light (EBL). We also want to determine physical quantities 

■ through the modeling of the object's broad-band spectral energy distribution (SED). 
^ ' Methods. VHE observations (~25 h live time) of the BL Lac 1ES 0347-121 (redshift z = 0. 188) were conducted with the High Energy Stereoscopic 
. • , System (H.E.S.S.) between August and December 2006. Contemporaneous X-ray and UV/optical observations from the SWIFT satellite are used 
^ • to interpret the SED of the source in terms of a synchrotron self Compton (SSC) model. 
JJTj ' Results. An excess of 327 events, corresponding to a statistical significance of 10.1 standard deviations, is detected from 1ES 0347-121. Its photon 
' spectrum, ranging from ~250 GeV to ~3 TeV, is well described by a power law with a photon index of T = 3. 10 ± 0.23 stat ± 0. 10 sys . The integral 
flux above 250 GeV corresponds to ~2% of the flux of the Crab Nebula above the same threshold. No VHE flux variability is detected within the 
data set. 

Conclusions. Constraints on the EBL density at optical to near-infrared wavelengths derived from the photon spectrum of 1ES 0347-121 are close 
to the strongest limits derived previously. The strong EBL limits confirm earlier findings, that the EBL density in the near-infrared is close to 
the lower limits from source counts. This implies that the universe is more transparent to VHE y-rays than previously believed. An SSC model 
provides a reasonable description of the contemporaneous SED. 
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1. Introduction 

1ES 0347-121 wa s discovered in the Einstein Slew Survey 
(Elvi s et al.l 1 19921) and later classified as a BL Lac object 
dSchachter et alj Il993b . Located at a redshift of z = 0.188 
dWoo et al.ll2005l) it harbors a sup er-massive blac k hole of mass 
log(M BH /M ) = 8.02 + 0.11 dWoo et all 12005b . The host is 
classified as an elliptical gal a xy with lumin o sity M r — -23.2 
dFalomo & Kotilainenl 1 19991) . IStecker et all (1 1996b used sim- 
ple physical considerations about the synchrotron and inverse- 
Send offprint requests to: M. Raue, e-mail: martin.raue@desy.de 
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Compton component of blazar spectra to calculate a prediction 
for the VHE y-ray flux above 0.3 TeV of 3.8 x 10~ 12 cnr 2 s" 1 
(0.03 Crab), which would be easily detectable with the cur- 
rent generation of VHE instruments. An upper limit on the in- 
tegral flux above an energy threshold of 1.46 TeV of 5.14 x 
10' 12 cm' 2 s-' (0.56 Crab) has been reported by the HEGRA 
collaboration (Aha ronian et alJl2004t) . considerably higher than 
the above mentioned prediction. 

For such a distant object, the observed VHE spectrum is ex- 
pected to be strongly affected by ab sorption of the y-rays due to 
the EBL dGould & Schredenll967l) . As a result, the VHE spec- 
trum measured here by H.E.S.S. provides constraints on the EBL 
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density (see, e.g., Aharonian et al.l l2006at hereafter Aha06) at 
near-infrared wavelengths. As 1ES 0347-121 is one of the most 
distant objects, for which a VHE spectrum is measured, these 
EBL constraints are potentially very strong. The broad-band 
SED of BL Lac objects, typically a double-p eaked shape, can 
often be explained by SSC models (see, e.g., [Aharonia n et al.l 
2005). However, the accuracy of such modeling is highly depen- 
dent on the observations being conte mporaneous due to the ex - 
treme flux-variability of BL Lacs (e.g. Krawczvn ski et al.1 l2000). 
The simultaneous UV/optical, X-ray, and VHE measurements 
presented here enable the first SSC modeling of the emission 
from this object including data in the higher-energy peak. 

2. H.E.S.S. observation and results 
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Fig. 1. The number of excess events versus squared angular dis- 
tance 9 2 to the nominal source position of 1ES 0347-121. The 
solid line shows the expectation for a point source from Monte 
Carlo simulations with the same pointing offset and zenith angle 
distributions as the observations. The vertical dashed line shows 
the cut on the distance from the source used for signal extraction 
(9 2 < 0.015 deg 2 ). 
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Fig.2. The differential photon spectrum of 1ES0347-121. The 
line is the best fit of a power law to the data (r = 3.10 + 0.23 sta t± 
0.10™). 



Fig. 3. The average nightly flux above 250 GeV from 1ES 0347- 
121. The dashed line shows the fit of a constant to the lightcurve 
(<£(£ > 250 GeV) = (3.32+0.44 stat )xl0- 12 cirT 2 s-'^/d.o.f. = 
30.4/28 = 1.1). 



The H. E.S.S. array o f four imaging atmospheric-Cherenkov 
telescopes (Hinton 2004) is used to search for VHE y-ray emis- 
sion from objects such as BL Lacs. The H.E.S.S. observations of 
1ES 0347-121 were performed between August and December 
2006. A total of 25 .4 h live time of good-quality data was 
recorded at zenith angles ranging from 12° to 40°. The mean 
zenith angle is ~19°, for which the energy threshold of this 
analysis is ~250 GeV. These data are analyze d with a standard 
Hillas-type analysis dAharonian et al.l l2006bl) . The data were 
recorded with a pointing offset of 0.5° relative to the nominal 
source position to allow a simultaneous estimation of the back- 
ground using e vents from the same fi eld of view {reflected back- 
ground from Ahar onian et al.ll2~006bl) . 

An excess of 327y-ray candidates is measured (1167 sig- 
nal events, 9241 background events, background normalization: 
0.0909) corresponding to a statistical significance of 10.1 stan- 
dard deviations (following Eq. 17 of iLi & Mal lT983). Figure Q] 
shows that the extension of the excess is compatible with that 
expected from a point-source. The fit position of the excess 
(a/2000 = 3 h 49 m 23.0 s ± 1.4^ t ± 1.3* yst , 6 jmo = -11°58'38" ± 
33" tet ± 20" yst ) coincides with the location of 1ES 0347-121 

(a/2000 = 3 h 49 m 23.2 s , 6 nm = -1 1°59'27.0"; ISchachter etafl 
1993). The angular distance between the fit position of the VHE 
excess and the source position is 47". The differential photon 
spectrum of the source is shown in Fig. [2fl A fit of a power- 
law function dN/dE = <f>o{E/l TeV)~ r to these data results in 
a statistically good description (x 2 /d.o.f. = 3.5/5) with normal- 
ization <D = (4.52+0.85 stat +0.90 sys )xl0- 13 cirrV 1 TeV" 1 and 
photonindexT = 3.10±0.23 sta t±0.10 S y S . The integral flux above 
250 GeV taken from the spectral fit is 0(£ > 250 GeV) = (3.9 + 
1-lstat) x 10~ 12 cirT 2 s _1 , which corresponds to ~2% of the flux 
of the Crab Nebula above the same threshold (Ahar onian et al.l 
2006b). As can be seen from Fig. [3] no significant variability is 
detected on time-scales of days or months. 



1 An electronic version of the HESS spectrum is avail- 
able from the publications section of the HESS website: 
http: //www.mpi-hd.mpg.c^7nfm /HESS/HESS.htmll 
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Fig. 4. The spectral energy distribution of 
1ES 0347-121. Shown are the H.E.S.S mea- 
sured spectrum (red markers), power-law 
fits to the intrinsic spectra corrected for the 
attenuatio n by different EBL models (solid 
red line - Primack et al.l 120051; dashed blue 
fine - P0.45 shape from [Aharonia n et al.l 
2006a), and contemporary SWIFT XRT and 
UVOT (black circles) and ATOM measure- 
ments (black triangles) (open black markers: 
contribution from the host galaxy sub- 
tracted) . Grey mar ke rs and li nes are archival 
(Eh 
JPe 
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A simple synchrotron 
self Compton model (see text) is also shown 
(solid curve). The model parameters are: 
6 = 25, B = 0.035 G and R = 3.2xl0 16 cm. 
The electron distribution is described by a 
broken power law with K = 1.26xl0 4 cm -3 , 

Trnin = 1X10 3 , y max = 3X10 6 , Threat = 

3.55xl0 5 , n, = 2.35, and n 2 = 3.35. The inlay 
shows a zoom on the intrinsic VHE spectra. 



3. SWIFT and ATOM observations and results 

SWIFT dGehrels et all 120041) observations of 1ES 0347-121 
were performed on October 3, 2006. A total of 3.2ksec of 
screened data in photon-counting mode are analyzed. For the 
analysis of SWIFT XRT and UVOT data the HEASOFT 6.2 
package with Xspecll and the latest calibration (XRT: 2007- 
03-30, UVOT: 2006-11-16) was used. The XRT data has been 
reprocessed with the updated calibration tables including new 
bad pixel tables. In order to avoid pile-up effects in the energy 
spectrum, the core pixels in the image were not u sed and the ex- 
tractio n region was chosen to be an annulus (e.g. Giom mi et all 
2006) with an inner radius of 5 pixels and an outer radius of 
30 pixels (1 pixel = 2.36")- The background spectrum was ac- 
cumulated in a wider annulus (inner radius of 45 pixels, outer 
radius of 90 pixels). The auxiliary response file was created 
with the standard tool xrtmkarf including the point-spread- 
function correction. A power-law model fit to the data (re- 
sponse matrix swxpc0to 12 J>00 10 10 lv008.rmf; photoelectric 
absorption fixed to the Galactic value of h# = 3.6 x 10 20 cm -2 ) 
between 0.3 and 8keV yields a good ;f 2 /d.o.f. = 105.3/105 = 
1.00 with r = 1.99 + 0.06 (error at 90% c.l.). The unabsorbed 
integrated energy flux between 2 and lOkeV from the model is 
/2-iokeV = (2.800 + 0.003) x 10- 11 erg cirT 2 s _1 . 

The Ultra Violet/Optical Telescope (UVOT) observations 
were made using 6 filter settings. Since none of the light curves 
from the different measurements indicate variability, the data 
sets for the individual filter wheel settings were added together. 
The aperture was chosen to be of 12 pixels radius for the op- 
tical filters while it is increased to 24 pixels for the UV filters 
(1 pixel = 0.48")- For the UV filter images, a 6-pixel radius re- 
gion centered on a faint stellar source 12" north of 1ES 0347-121 
falling into the wider-24 pixel aperture was excluded. The aper- 
ture photometry includes a pile-up correction and the flux has 
been calculated from the count rates by using the zero point val- 
ues quoted in the calibration notes of the SWIFT UVOT team. 
Finally, the observed flux was corrected for galactic absorption 



using a reddening of E(B — V) = 0.047 which was then scaled 
to obtain for the SWIFT filters following the recipe given by 
iGiommiet ail (120061) . 

The Automatic Telescope for Optical Monitoring (ATOM; 
iHauser etal .1120041) on the HESS site monitored 1ES 0347-121 
during the HESS observing period in November 2006. The mean 
B-band flux measured is 3.8 x 10 _12 ergcm _2 s _1 and the mean 
R-bandflux is 3.7xlO -12 ergcrrT 2 s _I (aperture radius: 4"). Both 
fluxes are constant, as all measurements are within 10% and 4%, 
respectively, of the mean values. 

The resulting SWIFT X-ray and ultraviolet-to-optical 
(UV/O) and ATOM data are shown in Fig. g] The X-ray and 
UV/optical fluxes observed are the highest measured from this 
source today, indicating a period of enhanced activity during the 
H.E.S.S. measurements. 

To correct for the contribution from the host galaxy in 
the UV /optical-bands the method presented in lAharonian et al] 
(2007) is followed. The host galaxy flux in the R-band of 
m r - 17.26 and t he half-light radius r e - 1.25" are taken from 
lUrrv et all (120001) . Using a de Vaucouleurs profile, the flux from 
the host galaxy falling in the signal aperture is estimated to be 
-80% for ATOM and -90% for UVOT. The host galaxy flux 
in the V, B and U-band are estim ated using the elliptic al galaxy 
spectral template at z — 0.2 from lFukugita et all (I 19951) . For the 
UVOT data host galaxy contributions of ~35% (V-band), ~1 1% 
(B-band) and -4% (U-band) are derived. For the ATOM data 
the host galaxy contributions are ~34% (R-band) and ~8% (B- 
band). The UV/optical measurements corrected for the contri- 
bution by the host galaxy are shown as open black markers in 

Rg. a 



4. Discussion 

The diffuse extragalactic photon field in the ultraviolet to far- 
infrared wavelength region (EBL) consists of the photons emit- 
ted by stars and absorbed and re-emitted by dust redshifted 
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Fig. 5. The spectral energy distribution of the extragalactic back- 
ground light (EBL). The black line is the limit (P0.61) on the 
EBL density derived from the measured spectru m of 1ES 0347 



121 u sing the same technique as described in Aharonian et al 
2006a). The dashed curve is the limit derived in lAharonian et al. 
2006a) from the HESS spectrum of 1ES 1101-232. Th e data 



compilation (grey markers) is from Mazinj&Raue] d2007l) . 



and integrated ov er time (see, e.g., [Hauser & Dwek 200l| and 
Kashlinsky 2005 for recent reviews). Its spectrum carries cos- 
mological information about galaxy and star formation history. 

Distant sources of VHE y-rays can probe the EBL density 
dStecker et alj|1992l) . VHE photons passing through the EBL 
are attenuated via pair p roduction: yvHE + Tebl — * e + + e~ 
dGould & Schrederlll967l) . Since this process is energy depen- 
dent, the VHE y-ray spectra measured on Earth carry an attenu- 
ation imprint from the EBL. With reasonable assumptions about 
the emission physics of the source, limits on the EBL density 
can be derived. 

Using the relatively hard energy spectra measured from the 
recently discovered VHE BL Lacs 1ES 1 101-232 (z = 0.186) and 
H 2356-309 (z = 0. 1 65), strong constraints on the EBL density in 
the optical to near-infrared (~0.8-4yum) were derived (Aha06). 
Following the exact same methodology as described in Aha06 
(scaling of a reference EBL shape until the intrinsic spectrum 
reaches a maximum hardness) a limit on the EBL density is de- 
rived from the observed VHE spectrum of 1ES 0347-121. For a 
scaling factor of 0.51 (i.e. a P0.51 shape), the fit of a power-law 
function to the intrinsic spectrum results in a photon index of 
r = 1.5, the minimum allowed value. Taking into account evo- 
lution of the EBL, again using the same arguments as Aha06, a 
limiting shape of P0.61 is derived. The limit is shown in Fig. [5] 
in comparison to the limit derived by Aha06 for 1ES 1101-232 
(P0.55). The limit derived here is only slightly less constrain- 
ing, which is a result of the softer observed spectrum and the 
smaller energy range towards low energies. The claimed excess 
of the EBL in the near-infrared abov e the value derived from 
source counts (Mat sumoto et alj|2005l) . often attributed to a pos- 
sible contribution of the first stars to the EBL, is again excluded. 

Given the large distance of the source, the measured VHE 
spectrum is severely altered by the EBL attenuation. Therefore, 
the VHE spectrum has to be corrected for this attenuation prior 
to any modeling of the intrinsic SED. Since the exact shape of 
the EBL is not known, the intrinsic spectr um is calculated for 
two different EBL models: the model from Prim ack et alj J2005) 



(PRIM), which has a low EBL density, and the upper limit shape 
P0.45 from Aha06 (AH). The intrinsic VHE spectra (inlay in 
Fig.[4]i are each well fit by power-law functions, which are shown 
in Fig.|4]above the measured VHE photon spectrum. The photon 
indices of the power-law functions are F = 2.10 + 0.21 for the 
PRIM and F = 1 .69 ± 0.22 for the AH shape. 

A simple homog eneous one-zone s ynchr otron self Compton 
(SSC) model from lKrawczvnskiet al. (2004) is used to describe 
the SED from the contemporaneous UVOT, X-ray and VHE 
data. For the modeling, the data are not strongly constrain- 
ing. Parameter sets describing the overall spectral shape with 
a standard shock-accelerated particle distribution can be found 
with a good statistical compatibility between the host-galaxy- 
substracted UVOT data, the X-ray data, both intrinsic VHE spec- 
tra and the model. In Fig.|4]an SSC model, fit to the host-galaxy- 
substracted UVOT, X-ray, and the intrinsic VHE data resulting 
from the AH EBL shape, is shown ~ 1.2; P(x 2 ) ~ 0.25). 
The archival radio measurements are assumed to be produced by 
a different particle population and are not included in the fit. The 
model parameters, although not strongly constrained, are similar 
to the parameters used previously to model other BL Lacs (e.g. 
lAharonian et alj|2007l) . 

To summarize, the new distant TeV blazar 1ES 0347-121 is 
detected at energies >250GeV with the H.E.S.S. experiment. It 
is among the most distant TeV blazars discovered today. The 
relatively hard VHE spectrum confirms previous limits on the 
EBL density at near-infrared wavelengths. X-ray and UV/O data, 
taken during the H.E.S.S. observation period with the SWIFT 
satellite, indicate a high flux state of the source. The contempo- 
raneous X-ray, UV/O, and VHE data are not strongly constrain- 
ing and can be described by a simple SSC model. Further multi- 
wavelength observations, especially in X-rays, can improve the 
constraints on the model parameters. Many TeV blazars show 
large variations in the VHE flux, so monitoring the source, to 
measure a state of high VHE flux (flare), is desirable. 
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